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Abstract 

We have studied one-dimensional cellular automata with updating rules depend- 
ing stochastically on the difference of the heights of neighbouring cells. The prob- 
ability for toppling depends on a parameter A which goes to one with increasing 
slope, i.e. the dynamics can be varied continuously. We have investigated the 
scaling properties of the model using finite-size scaling analysis. A robust power- 
law behavior is observed for the probability density of the size of avalanches in a 
certain range of A values. The two exponents which determine the dependence of 
the probability density on time and system size both depend continuously on A, 
i.e. the model exhibits nonuniversal behavior. We also measured the roughness 
of the surface of the sandpile and here we obtained an universal behavior, i.e. a 
roughness exponent of about 1.75 for all values of A. For the temporal behavior 
of the mass a spectrum is obtained with an exponent </> close to 2 again for 
all A- values. 



1. INTRODUCTION 

Bak, Tang and Wiesenfeld |]J recently introduced the notion of self-organized criticality 
(SOC) as a framework to understand the dynamics of extended dissipative dynamical 
systems. The dynamics drive the system into a state where relaxation processes occur 
over all time and length scales. It was suggested [l[] that sandpiles were a particularly 
clear example of a self-organized system. 

Theoretical sandpile models have been investigated by several authors Q-f|. They 
usually involve two kinds of steps. First a particle is added at ramdomly choosen sites, 
then the particles are eventually distributed to neighboring sites according to dynam- 
ical rules. Usually, these rules are deterministic: there is toppling of particles or not 



depending on the local occupation of sites. However, real sandpiles behave in a different 
fashion. In real sandpiles the particles have different shapes, various sizes and they can 
stick together. As a matter of fact, large and very unstable slopes can develop with 
slides of particles in an apparently stochastic fashion: a certain slope may or may not 
lead to a slide depending on a great number of internal degrees of freedom like the 
above mentioned different shapes and sizes of the particles, the friction between them, 
the local variations in packing etc. To incorporate some of these features into a the- 
oretical model we propose as a first step a model with dynamical rules which are not 
deterministic: toppling of particles occurs with a probability function which is small 
if the slope to neighboring sites is small and which increases to one with increasing 
slope. In the present paper we report results of extensive simulations of this model. We 
restrict ourselves to the one-dimensional case which shows already an interesting and 
nontrivial behavior. 

2. THE MODEL AND SIMULATIONS 

The model we consider is an one- dimensional "nonlocal unlimited" sandpile model 
p| with the following rules. We assume integer heights hi at lattice sites % = 1, 2, 3, L, 
where L is the size of the system. The local slope Oi at the site % is defined as the 
height difference between two neighbouring cells, cr^ = hi — h i+1 . The dynamics obey 
the following rules: If the slope cr, exceeds a threshold value a c an avalanche can be 
generated. With a certain probability p(<Ji) grains drop from site i to the neighbouring 
cell with 

h{ > hi Tii 

h i+j -> hi+j + 1 for j = 1,2,..., m (1) 

The number rij of grains which topple grows with increasing <7i, rii = <Ji — N (unlimited 
model). Since there is toppling to rii right neighbours the model is called nonlocal. 
Note that from Eq.(l) it follows that all slopes in the system are non-negative. 
The simulation starts by adding particles at the top until the stability condition o\ < a c 
is violated. Then, with probability p(eri), n\ grains flow to the right neighbour according 
to Eq.(l) and an avalanche develops. If no toppling takes place another grain is added 
to site 1. If toppling takes place the next site will be visited and so on. An avalanche 
stops at site i if this site doesn't topple or if its slope to site i + 1 is too small. Thus an 
avalanche only starts at the top. The length of an avalanche is defined as the number 
of sliding sites. At the right edge of the system grains can leave the pile, i.e. hi is set 
equal to zero for % > L. 

The most important new ingredients in our model is the probability function p(crj) 
which determines whether a certain site i topples or not. This function should obey the 
following conditions: 

p(a < cr c ) = 

p(a ) > p(a) for a > a (2) 
p(a) — > 1 for a — > oo 



The first condition states that there is only toppling if a threshold value is reached, the 
second conditions means that there is increasing probability for toppling if the slopes 
increase and the third equation means that toppling always takes place for very large 
slopes, i.e. the heights hi are finite for any system of finite length. In our simulation 
we use 

p(<Ti) = 1 - e-^- 1 ) (3) 

A is called the toppling parameter. Thus there are three parameters in our model: a c 
which determines the stability condition, N which determines the number of sliding 
grains in one event and most importantly the toppling parameter A with which we can 
tune the dynamics continuously. In our simulations we choose a c = 2, N = 1 and study 
in particular the dependence of the model on A. Note that for A —>■ oo one gets a trivial 
model: erj = 1 for all i and any avalanche reaches the edge of the pile. 
The most important quantity to study is the probability density for the lengths of the 
avalanches for which a finite-size scaling analysis is performed. As usual we assume for 
this probability density P(t, L) a scaling form 

P(t,L) = L-?g(L-»t) (4) 

where t is the number of toppling sites. For t <C L we expect a power-law behavior, 
P(t, L) ~ t~ K independent of L. Thus (3,u and n must obey the following scaling relation 

@ 

(3 = uk (5) 

Another important quantity is the roughness exponent ( which a is measure for the 
fluctuations about the average height profile and which reflects the critical nature of 
the steady state. It is defined by ||: 

^E<(^-<^>) 2 > (6) 

i=i 

Finally we measure the total mass of the pile 

M(t)=J>(t) (7) 
and subsequently calculated its Fourier spectrum |M(/)| 2 to determine the exponent 

4>. 




3. RESULTS 

The simulations were performed for lattice sizes of L = 50, 100, 200, 500, 1000. Before 
performing a measurement a large number of particles (arround 10 8 ) was added start- 
ing from an empty lattice to reach the critical state. Note that all measured values 
are independent of the initial conditions, i.e. the critical state is an attractor of the 
dynamics. 

First we measured the probability density of the length of the avalanches P(t, L) for 
a lattice size L = 100 and examined how P(t, L) depends on the toppling parameter 
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Figure 1: Finite-size scaling fit for A = In 5 and five different lattice sizes L = 
50, 100, 200, 500, 1000. The solid line corresponds to a power-law with an exponent k = 2.453. 
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Figure 2: Scaling exponent (3 and power-law exponent n as a function of the toppling param- 
eter A. Solid circles correspond to v~ l (3. Eq. (5) is fulfilled within numerical accuracy. 



A. A power-law behavior is observed in the range 1 < A < 3. Outside this range we 
obtained on a log-log plot curves which deviate visible from a power-law behavior. We 
have not investigated this problem any further but restricted ourselves to the A-values 
where power-law behavior is clearly observed. Then the probability density P(t, L) 
was measured for various values of A and lattice sizes. As can be seen from Fig.l 
finite-size scaling works extremely well. For different values of A we get v as 1 but the 
scaling exponent (3 changes strongly [Fig. 2]. Thus we expect a varying exponent k. 
k is calculated using regression analysis for lattice sizes L = 100,200,500, 1000. The 
averaged value of k is shown in Fig. 2. Note that the change of the exponent can not 
be explained by statistical errors since they are far too small. The relation between the 
exponents, Eq.(5), is always fulfilled. 

To determine the roughness exponents ( we limit our simulation to lattice sizes L < 500. 
After reaching the critical state we add around 10 8 particles to calculate £ for different 
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Figure 3: Roughness £ as a function of lattice size for different values of the toppling parameter 
A. The solid line corresponds to an power-law £ ~ Z> with £ = 1.75. 
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Figure 4: Power spectrum of the time evolution of the total mass of the pile for A = In 5. 
The spectrum is an average over 40 measurements with 2048 events. Solid line corresponds 

to \MU)\i~r 2M9 - 



values of A. Fig.3 shows the results. One can see that for L > 100 the roughness £ scales 
with the system size as and ( is independent of the toppling parameter A. However, 
one should note that for a very accurate determination of ( one needs simulations with 
system sizes larger than 500 since the data points shown in Fig.3 are still not on a 
straight line. But those measurements would demand more than 10 9 events and have 
not yet been done. Increasing the lattice size the roughness grows very fast, because 
C > 1. The average height profile scales as < hi >~ L ai with a± ~ 1.8 for all i and A 
Q The exponents a± and C should obey the relation ot\ > ( which is fulfilled in our 
model. It is very interesting to see whether the roughness exponent ( is equal to the 
scaling exponent of the averaged height a±. Finally we have calculated the time depen- 

1 Notice that the average height profile is not a linaer function of i, but obeys the scaling relation 
< hi >= L ai f(L- a H) with f = a + a ie - a2%2 §. 



Table 1: The critical exponents for different values of the toppling parameter A. 
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dence of the total mass of the pile (Eq.7). We used a lattice size L = 50 and measured 
the mass M{t) for 2048 events. Fig. 4 shows the average over 40 of the corresponding 
fourier-spectra. The power spectrum |M(/)| 2 scales as with ps 2 for all A. 

4. CONCLUSIONS 

In conclusion we have simulated an one-dimensional cellular automaton with nonlocal, 
unlimited and stochastic dynamics which exhibits SOC behavior in a certain range of A- 
values. All trademarks of SOC are observed: finite-size scaling and a power-law behavior 
of the avalanche probability as well as behavior for the power spectrum. However, 
our model displays nonuniversality: the exponents (3 and k change continuously with A. 
On the other hand the exponents which describe the fluctuations, <p and (, are universal, 
i.e. independent of A. Kertesz et al |J have shown that a great class of sandpile 
automata yield f~ 2 behavior. Thus it is understandable why <fi is constant. In contrast 
reasons why the roughness exponent is independent of the toppling parameter are not 
known to us. Generalization of this model to higher dimensions and an investigation 
of the behavior for very small and very large values of A where power-law behavior 
apparently breaks down is under way. 
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